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ABSTRACT 



We present a geometrical methodology to interpret the periodical light curves 
of Soft Gamma Repeaters based on the magnetar model and the numerical arith- 
metic of the three-dimensional magnetosphere model for the young pulsars. The 
hot plasma released by the star quake is trapped in the magnetosphere and 
photons are emitted tangent to the local magnetic field lines. The variety of 
radiation morphologies in the burst tails and the persistent stages could be well 
explained by the trapped fireballs on different sites inside the closed field lines. 
Furthermore, our numerical results suggests that the pulse profile evolution of 
SGR 1806-20 during the 27 December 2004 giant flare is due to a lateral drift of 
the emitting region in the magnetosphere. 

Subject headings: gamma rays: bursts - stars: magnetic fields - stars: neutron - 
X-ray: stars 



Introduction 



Soft Gamma Rep eaters (SGRs) seemed weird since the first discovery on 1979 March 
5 (IMazets et al. Ill979h for their mysterious characteristics such as the large energy release, 
the repetitive emission of bursts in hard X-rays or soft gamma-rays bands, and the pulsed 
periodical emissions after the bursts and in the quiescent stages whose m orphologies are both 
energy-dependent and time-dependent (see lWoods fc Thompson I (120041 ) for a recent review). 
So far, the catalogu^ has four SGRs confirmed plus one candidate. SGRs are found to be 
associated with young (~ 10'^ yr) supernova remnants (SNRs), and their spin periods are 
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5 ~ 8s and at a large spin down rate about 10 ^^s s ^, which give an inferred ultra-strong 
dipolar magnetic field at the order of 10^^ ~ lO^^G. 

A variety of models were proposed to understand the physics of SGRs, and the magnetar 
model i s now widely accepted to asc r ibe SG Rs as neutron stars with magnetic field of 10^^ ~ 



lO^^G (IThompson fc Duncan Ill995l . Il996l ). Unlike most of the pulsars in the neutron star 
family powered by their spin-down, the high- luminosity bursts and the persistent X-ray 
or soft 7-ray pulsations of magnetars come from the decay of their ultra-strong magnetic 
fields. The spectrum of the persistent X-ray emission could be fitted by a superposition of a 
blackbody component and a power law, which suggests that besides the radiation from the 
neutron star surface, there is another component coming from the magnetosphere. 

The high-luminous burst has now been successfully interpreted by the dissipation of 
magnetic energy. However, the persistent long-period pulsations in the quiescent X-ray 
emission o f SGRs are still not well understood for their complicated and astonishing mor- 



phologies. [Thompson. Lyutikov fc Kulkarni I (120021 ) assumed that the angular pattern of the 



X-ray fiux are modified by the resonant cyclotron scattering at a distance about 50-100 km 
above the neutron star, and multi-pulses are generated by some effects associated with the 
twisted magnetosphere (e.g., the optical depth is thin near two magnetic poles, and thick at 
the magnetic equator). However, their numerical results from the Monte-Carlo simulation 
has a bias in favor of the orthogonal dipole (i.e., the magnetic axis is perpendicular to the 
spin axis), and the viewing angel h as also to be nearly 90° from the spin axis, for the photons 



escaping from optically thin poles (IFernandez fc Thompson 1120071 ). In addition, the separa- 



tion of the strong peaks in their calculations is always tt (half phase cycle), which differs from 
the observation. Thompson (private communication, 2005) agreed that the change in the 
persistent pulse profile refiects a re-distribution of persistent currents on closed field lines. 
However, he ascribed such re-distribution as the gradual change in the magnetic twist at a 
distance of 30-100 km, while our understanding for that is due to the gradual migration of 
the emitting regions (crustal platelet motion or hot plasma drifting in the magnetosphere, 
we will discuss the possibility of each candidate in 

In this paper, we present some simulated pulsed profiles by assuming the radiation region 
located in the closed field lines, and make attempts to simulate the radiation morphology 
evolution in one particular event, i.e., SGR 1806-20 burst on 27 December 2004. We first have 
a general review on the timing properties of SGRs in ^ In ^3], we present the motivation 
to reproduce the light curves in the closed field lines and the resultant profiles. In ^ we 
applied the geometrical model to the well-known December 2004 burst of SGR 1806-20 and 
calculate the radiation morphology by a three-dimensional magnetosphere simulation. We 
also try to explain the change in the persistent pulse profile. A brief discussion is given in 
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Periodical Emissions from SGRs 



The periodical emissions of the four confirmed SGRs have been detected in the decay 
of burst and the quiescent or persistent stages, with the periods in the 5-8 seconds range 
which are the spin periods of the neutron stars. The pulse profiles show many interesting 
and even surprising morphologies, some of which are totally different from the canonical 
radio or high energy pulsars. The characteristics of light curves of the X-ray and 7-ray 
pulsars such as number of peaks, the peak separation and the relative amplitudes of the 
peaks are unchangeable, while those of the SGRs are time dependent. Here, we list several 
main features of the persistent emissions from SGRs: 

1. Multi-peaked morphology: The most dramatic example is the pulse profile change 
of SGR 1900-1-14 after the outburst on 1998 August 27. A four-peaked repetitive pattern 
of the X-ray lig ht curve was detected by both Ulysses and BeppoSAX half a minute after 
the burst onset (IHurley et al. Ill999l : iFeroci et al. Il200ll ). and these peaks were found to be 
evenly spaced at 1.0 s intervals on the 5.16 s rotation period. A nd several minutes later, thi s 
multi-peaked profile evolved as a simple sinusoidal morphology (IWoods fc Thompson 112004 ). 
and this evolution in pulse profile lasted for a couple of years. Such change in morphology is 



also found in the RXTE PGA archive of SGR 1806-20 between 1996 and 2005 ( Woods et al. 



20071 ). During the 10- year observation, the 2-10 keV pulse shape evolved from one broad 
peak pattern to a three-peaked one in 2003, and then the sinusoida l shape again until t he 
multi-peaked profile after the burst on December 2004 (see Fig. 3 in lWoods et al. I (120071 )). 



2. Relative magnitudes of peaks: Besides the evolution of the number of peaks, the 



relativ e magnitudes of the peaks in one phase cycle also change with time. iPalmer et al 



( 120051 ) showed such pulse profile evolution during the giant flare of 27 December 2004. The 
folded light curves in different time intervals from 30 to 265 seconds following the main spike 
indicate the growth of the second peak, whose intensity related to the primary peak increases 
from the DG level to nearly equal in height. In the other words, we can say that the primary 
peak fades until the same magnitude as the secondary one. At the late stage of the decay 
of the giant flare, the relative magnitude of the third peak 0.2 in phase prior to the primary 
one starts to grow up. 

3. Energy- dependent profiles: The evolution of p ulse profiles of SGRs is not only time- 
dependent, but also changes in different energy bands. IWoods et al. I (120071 ) investigated the 
energy dependence of the SGR 1806-20 pulse profiles in three energy bands between 2 to 40 
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keV. Six months before the flare, there was only one broad peak in the pulse profile below 
15 keV, and it showed two clear peaks in the 15-40 keV band. After the giant flare, the 
light curve became more complicated, showing m ultiple peaks in all e nergy intervals, and 



the peaks were inconsistent in phase (see Fig. 4 in lWoods et al. I (120071 )) 



The totally different phenomena require the totally different physics in this small (may 
be not small) group of neutron stars, compared with the canonical pulsars. The thermal 
spectrum component suggests that the radiation partly comes from the hot spots on the 
stellar surface. In addition, the phase inconsistence of the pulse profile indicates that it may 
not be localized in a particular region, and the non-predictable star quakes can produce the 
randomly- localized emitting regions during every burst event. This is the main motivation 
for us to build up the model of the alterable pulse profiles in the next section. 



Radiations from the Closed Field lines 



3.1. Why closed regions? 

The theoretical models for radiation from high energy pulsars (e.g. Crab and Vela) re- 



polar gap model ( 


Harding 


1981 




Dauffhertv & Harding 


199e 


), or the outer gap model 


(Chene. Ho & Ruderman Ill986: 


Chiang & Romani 


1994: 


lomani & Yadie 


aroelu 1995: Chene. 


2000 


), or the modified outer gap model ( 


Dvks & Rudak 


2003; 


Jia et al. 


2007 


). However, 



we cannot apply them to magnetars to explain their persistent X-ray emissions. We assume 
that the pulsed emissions of SGRs in the decay or afterglow of the bursts are originated in 
the closed regions rather than in the open regions. We have three main reasons for confining 
the radiation regions to the closed zones in the neutron star magnetosphere: 

1. The volume occupied by the open field lines is much smaller compared to that by the 
closed field lines. Take SGR 1806-20 for example, the spin period is P = 7.56 s, which means 
the radius of the light cylinder reaches Rl = f~ = 3.6 x 10^° cm, where c is the speed of light. 
Thus, the radius of the polar cap is Rp = R^^^/R^/Rl = 5.3 x 10^ cm, corresponding to an 
angular size of 6p = 0.3°, much smaller than the characteristic size of the crustal platelet, 
where R^, = 10^ cm is the stellar radius. So, the polar cap area is only 10~^ of the neutron 
star surface, and there is no reason why the crustal platelet, where the magnetic energy is 
released, should be located in the open area. 

2. The open field lines reach the light cylinder, where the co-rotating speed approaches 
the spee d of light, and the relativistic effects play a significant role in the radiation mor- 
phology (IRomani fc Yadigaroglu Ill995l ). Such effects lead to the sharp and narrow peaks 



- 5 - 



of the light curves, and both peaks are produced by one single pole. However, those effects 
become less important when it is applied to the closed field lines, for they are much closer 
to the stellar surface than the open lines. As shown in Fig{T|, we find that the farthest 
distance where the closed field lines can reach drops much quickly when their footprints are 
displaced further away from the polar cap. For the closed line on the plane (on which both 
the magnetic axis and spin axis lie) with layer parameter oi = 5, it only reaches a distance 
less than 0.05/2/, away from the stellar surface, and O.OIRl for oi = 10 (the definition of Oi 
is given in §3.21) . Thus, the double-peaked profiles are not necessary for closed field hues, 
and the broad peak could be a general product instead of the sharp one. Other features, like 
the peak separation and the multi-peaks, could also be explained by the closed field lines 
(details to be discussed in the next section). 

3. The emitting region has not to be localized on some particular sites on neutron star 
surface or inside the magnetospher e while in the open field lines, the accelerating gap is 



restri cted to the null charge surface (iGoldreich fc Julian Ill969l : ICheng. Ruderman fc Zhang 



2000l ). This freedom in loci ensures the variety of the pulse profiles during different outbursts, 
for we can argue that time-dependent light curves result from different bunches of magnetic 
field lines where the plasma is trapped. 



3.2. Strategies for numerical simulation 



Since we argue that the periodic pulsed radiation of magnetars in the quiescent state 
is released from the close d field line regions rather than the open one, we could apply the 
3D magnetosphere model (jCheng. Ruderman fc Zhang 1 120001 ) to simulate the pulse profiles. 
The boundary of these two different regions is defined as a bunch of so called the last closed or 
the first open magnetic field lines, which are tangential to the light cylinder. The hot plasma 
is trapped in the closed region and oscillates along the closed field hues, and photons are 
assumed to be emitted outwardly along the tangential direction of the magnetic field lines, 
which makes the pulsed radiation much different from the one generated in the open area. In 
order to be consistent with the pulsar magnetosphere calculation in the high-energy pulsar 
models, we adopted the same definition for the coordinates of fo otprints of the magnetic 
field lines and layer parameters ljCheng. Ruderman fc Zhang Il2000l ). The shape of the polar 
cap could be determined by the footprints of the last closed field lines anchored in the stellar 
surface, and we can label the coordinates of these footprints as (xq, yo, Zq). Then we are able 
to define another set of footprints of magnetic field lines (xq, ?/q, Zq) by multiplying a factor 
ai called layer parameter: x'q 



aiuo, and z'q = ^/l- (a;[f + y'^), where ai > 1 



indicates the closed region, and ai < 1 represents the open ones. 
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Since the magnetosphere is co-rotating with the neutron star, aberration effect occurs 
along the line of sight. Thus, we have the aberrated emission direction (in the observer's 
frame) u' = {u^,Uq,u^ in the expression of the direction u = {ur,U0,u^) in the co-rotating 
frame and the rotational speed f3 = \r x fl\/c: 

~ (1 + ^u^c) 
(1 + f3u^c) 



Ur 



(1 + (3u^c)' 

Another effect taken into account is the time of flight, which differs a lot for the photons 
originated at different sites inside the magnetosphere. This effect can lead to the phase 
difference of the arrival photons comparable to the light curve period. Combing these two 
effects, and choosing th e rotational axis as the 2;-axis, we obtain the phase angle $ and the 
polar angle ( given by (jYadigaroglu 1119971 ) 

cos( = u^/u , 

$ = -0\-r-^7i?i, (2) 

where 0\ = aiccos{u^/u^y) (here (x, y, z) is the cartesian coordinate system) is the azimuthal 
angle in the observer's frame. Choosing Q,-^ plane to be the x-z plane, u'^y is the length of 
the projection of ii' on the x-y plane. 



3.3. Simulated light curve profiles 

In the following, we adopt in this paper the magnetic coordinates (0*, (j)^) to describe the 
loci of radiation regions in the magnetosphere. The polar angle 9^ is defined as the angle with 
respect to the magnetic axis, instead of the rotational axis. The emitting region is assumed to 
be in the shape of a band along the azimuthal direction, with a relatively smaller longitudinal 
thickness compared with its azimuthal width. As the neutron star is rapidly rotating, we 
set the phase of the Vt-fx plane defined by the rotational axis and magnetic axis as 0* = 0°. 
Therefore, the transverse extension of the emitting region along the azimuthal direction could 
be expressed as which is treated as a parameter in our numerical simulations in this 
paper. However, there is no particular parameter for a quantitative longitudinal thickness 

in our simulation, and we represent it by the layer parameter ai, instead. 

Combining those factors mentioned above, we show in Figl2] and |3] two typical light 
curves generated in the closed field line zone. Figl2] is a single sinusoidal pattern, and Figj3] 
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is a double-peaked morphology without any off-pulse phase. The spin period of SGR 1806-20 
is applied in the calculation, and the inclination angle and viewing angle w.r.t. the spin axis 
are chosen to be 3 0° and 50°. Bo t h puls e profiles are commonly detected in the SGRs timing 



observations (e.g. IWoods et al. I (120071 )). The azimuthal widths of the emission regions for 



these two cases are A0* = 180°, and the other parameters for these two plots, like oi and 
the (6'*,0*) coordinates of the emitting regions, are given in TableiH In the upper panels 
of both figures, we show the emission projections onto the {(, $)-plane to illustrate the two 
kinds of light curves, respectively. We can find that the whole {(, $)-plane is fully filled by 
the emissions from closed fiel d lines, compared with those partiall y-filled by the outer gap in 
the open lines (e.g. Fig. 6-8 in lCheng. Ruderman &: Zhang I (2000)). The emission projection 
of Figl3]has a dense bundle in the phase range (0.5, 0.7), which produces the secondary peak 
of the double-peaked light curve. 

What makes such difference b etween these two kinds of light curves? As shown by 
Cheng. Ruderman fc Zhang I (120001 ). the neutron star rotation results in a non-uniform dis- 
tribution of the magnetic field lines in the magnetosphere, i.e., the magnetic field lines are 
swept back to accumulate around the Q-fi plane. In our calculation, the only different pa- 
rameter in both cases is the longitude of the emitting region. On the side around 0,, = 0°, 
the accumulated field lines produce the one-peak dominated light curve profile (e.g. the 
single sinusoidal pattern), while on the other side around 0* = 180°, the widely-separated 
field lines make the double-peaked profile possible. 



4. Application to SGR 1806-20 in December 2004 

4.1. Pulse profile evolution 

As illustrated in the standard magnetar model, the crust of the magnetar breaks when 
the magnetostatic equilibrium in the lower crust is no longer be sustained, and launches 
a hot fireball, which triggers the outburst of SGRs. The released energy comes from the 
reconnection of magnetic field lines in a crustal plate, which can be modelled as 

E^^l\ (3) 

where B is the magnetic field in the crust, and I is the size of the crustal plate. The energy 
released by SGR 1806-20 in the Dec 27 2004 burst is estimated to be as high as ~ lO^^erg, 
and by substituting the inferred magnetic field strength of order 10^^ G, we can estimate 
I ~ 10^ cm, which is about the thickness of the neutron star crust. A clump of electron- 
positron or electron-proton plasma is then ejected into the magnetosphere, and trapped by 
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the magnetic field lines anchored in such crustal platelet. The charged particles emitted 
by the hot plasma travel along the closed magnetic field lines, and radiate photons. In the 
following emission morphology simulation, we will assume the emissivity along the field lines 
is uniform. However, the emission region along the azimuthal and polar directions have 
a finite characteristic dimension correspondi ng to ~ lO^cm. We assumed that either the 



crustal motion driven by the neutron vortex (IRuderman Ill99ll ) or the lateral motion of the 



plasma across the field lines driven by the residual electric field, could lead to the change of 
the radiation morphologies. We will discuss which mechanism is more plausible later. 

We attempted to simulate the pulse profile evo lution during the giant fiare of 27 De- 
cember 2004 (e.g. Figure 2 in Palmer et al. (j2005 )). The most significant feature of the 



SGR 1806-20 pulsed radiation is the increased amplitude of the secondary peak related to 
the primary one. Our trials on the numerical simulation suggest that such change is due to 
the motion of emitting region. We present our results in FigJU which shows the effect of the 
azimuthal motion of the emitting region across the magnetic field lines with the typical layer 
parameter oi = 10, anchored in the crustal plate with the size about 1.5 x 10^ cm. FigHlalso 
shows some other features (e.g. the width of the peaks and the phase separation between two 
peaks) to be consistent with the observation. The inclination angle of the magnetic dipole 
and the viewing angle of the observer are chosen to be 30" and 50^, respectively. We assume 
that size of the plasma or the emitting region remains unchanged during the motion, and 
we describe the loci of the emitting region center in terms of the {6*c,4'*c), which are the 
magnetic coordinates of the middle point of the emitting region at the typical layer. Since 
the motion is along the azimuthal direction, 6^^ = ai x 6p = 3^ for all three panels in FigJH 
As the Q-fi plane is defined as 0* = 0°, the azimuthal angle 0* increases along the spin 
direction of the neutron star. The center of the radiation region was shifted from 0*c = —11*^ 
(panel a in FigHj) to 0=kc = 0° (panel b in FigJl]), and finally at 0=kc = 7.5° (panel c in FigJlj) 
All parameters applied in our model fitting are listed in TableO 



As shown by lArendt &: Eilek I (Il998l ). the polar cap of a rotating neutron star is asym- 
metrical and even probably discontinuous, which could make the radiation morphology com- 
plicated and asymmetrical. In our simulation, the primary peak at phase 0.9 of the light 
curve results from the majority of the magnetic field lines in which the plasma is trapped, 
and only a small portion of the field lines on the right edge of the plasma generate the 
secondary peak (at phase 1.2), where 'right' means the site whose 0* value is larger. As 
the plasma is drifted from left to right (along the direction 0=,, increases), more magnetic 
field lines at the right site are enrolled to radiate photons, which makes the secondary peak 
grow up. The radial distance of the local place where the arrival photons are generated is 
shown in FiglHl Since the magnetosphere co-rotates with the neutron star, we calculate the 
modelling velocity of the emitting region motion in (p^ direction. As shown in FigjTl the 
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speed of the drift is estimated to be ~ 10"^ cm s"^ according to the pulse profile evolution 
timing by observation. We also give the expected radiation from the two magnetic poles in 
panel c of Fig. HI which are indicated by two arrows. As the definition of the {Q, /i)-plane, 
the magnetic poles are located at 0* = and vr (0 and 0.5 in the phase cycle), respectively. 
However, the time of flight effect makes the radiation from magnetic poles has a tiny shift in 
phase, e.g. 0.1 and 0.6. At the early stage, the flux of the radiation from the trapped plasma 
is so strong that the emissions from the polar caps could not be resolved. However, as the 
intensity of the two main peaks fades, the two minor peaks might become discriminable. 

In general, the emitting region should not always move in one direction. Fig. [5] shows 
the pulse profile evolution when the motion of the radiation region is changed to be along 
the ^-direction. Panel c is the same as the one in Fig. HJ and the light curve of panel d is 
produced when the center of plasma is located at ai=20, with the same azimuthal position 
as that in panel c. We find that the morphology of the pulsed emission remains roughly 
unchanged, corresponds to the late stage (~ 170s after the burst) of tail in the giant flare of 
27 December 2004. We want to remark that the light curve can also remain unchanged when 
the motion of the radiation region stops. We cannot differentiate these two possibilities from 
the light curve evolution. 



4.2. Interpretation of numerical results 

In order to simulate the time evolution of the light curves of SGR 1806-20, we have 
assumed that the emitting region can migrate in the magnetosphere. In the following, we 
would like to discuss several possible movements induced in the crust and the magnetosphere, 
and their resultant speeds. 

The interaction between the flux tubes and vortex lines in the core of neutron star make 
them inter-pin to each other. The vortex lines will move out of the core due to the spin-down 
of the star, they drag the flux tubes with them. However, the flux tubes are anchored in 
the crust, consequently a large stress will be applied to the crust from the flux tubes (e.g. 
Ruderman 1991). When the crust breaks, the flux tubes will try to reduce their tension and 
drag the broken crust platelet to move. The maximum shear stress on the base of the crust 
from core magnetic flux tube motion which the crust could sustain is 

S ^ ^ 3 X 10^9 dyn cm'^, (4) 

Svr 

where Bcru ~ 10^^ G is the critical magnetic field inside the magnetic flux tube. However, 
the crust may already break before the magnetic stress reaches the maximum value. The 
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shear modulus yU of a 1 km thick crust could not be larger than 10^° dyn cm ^, and the crust 
may break when the stress reaches 



S: 



break 



(5) 



where / is a factor of order unity, and 6max is the yield strain under tension or compression, 
which is ~ 10~^ — 10~^ (c.f. Ruderman 1991). 

Another mechanism to drive the crust to move is the vortex creeping. As the rotation 
of the neutron star slows down, the flux tubes are driven outwardly by the the neutron 
vortices. The force acting on unit length of a flux tub e at the core-crustal interface is 
jchau. Cheng fc Ding I[l992 : IPing. Cheng fc Chau Ill993h 



fn 



(6) 



where $o = hc/2e ~ 2 x 10^''G/cm^ is the flux quantum, and the subscript c represents the 
values in the core. VLg is the rotation rate of the core superfluid, and can be approximated 
as that of the crust Vt^ or the spin rate of the star in our following estimation, uocr is the 
maximum angular velocity lag between Vt^ and ^7c. is the core magnetic field, and could 
be only about 10~^ of the surface value because the flux tubes are pushed out of the core 
due to spin-down (Ding, Cheng & Chau 1993). Thus, we can estimate the total magnitude 
of the force acting on the crust platelet with area A. The number of flux tubes anchored in 
this platelet is given by 



iV/ 



total ' 



A TlRlBr 



AB, 



Substituting /„ in Eqj6|, we have the total driving force on a platelet 

Fn = NffnRc = 2ApcRlfls^cr- 

We then make a dimension analysis to estimate the velocity 

<v>r^ VF/pA = {Rl^l,uj,rf/\ 



(7) 



(9) 



ujcr is given by various models as of order 10 ^ rad s ^ (IChau. Cheng &: Ding Ill992l ). there- 
fore, we have the velocity of the flux tube about 10^ cm s~^. 

However, when the flux tubes move, they will experience a drag force by the electron 
sea in the core. The drag force per unit length of a single flux tube is given by 



64 EfKc 



(10) 
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Here, Ue ~ 10 cm~ is the electron density, and Ef is the electron Fermi energy, which is 
about 100 MeV. The penetration length of a proton A is ~ 100 fm, and Vp the velocity of 
the flux tube. By equating the vortex acting force and drag force fy, we calculate the 
velocity as 



64 EfKcpcRc^s^c 
10"^^ cm s"\ 



(11) 



and iRuderman. Zhu &: Chen I (119981 ) have considered the collective motion of flux tubes and 



gave an estimation of 10 ^ cm s ^. If we equate the magnetic stress force SA and the drag 
force by 



BB, 



we have the velocity 



A = Uvp)NfR, 



64 BEfkc 



(12) 



(13) 



at about lO^'^ ~ 10"^ cm s ^. Therefore, we find the velocity of the crustal motion is too 
small to account for the emitting region drift. 

On the other hand, if we consider the drift velocity of the plasma driven by the electric 
field in the magnetosphere, we can write 



5E X B 



52 



-c ~ 



5E 



B{D) 



c, 



(14) 



where 5E is the residual electric field which drives the plasma to move, and B{D) is the local 
magnetic field at the radial distance D, where the emitting plasma is located. As shown in 
FigEl the average radial distance of the emitting region is about D ~ 0.002Rl = 0.7x 10^ cm. 
The length of the magnetic flux loop at ai = 10 is of order about C ~ O.OIRl referred to 
FigUl The residual electric field could be determined by 6E ^ AV/C, where AV^ is the 
electric potential drop. 

In closed field it is not clear how any substantial potential can survive because elec- 
trons and positrons can be created and trapped in the closed field line region. These 
electron/positron pairs can screened the electric field. In our simulation, we find that 
the emission region is characterized by ai=10, which is not far away from the open field 
line region. For a platelet with characteristic dimension of ~ lO^cm, we can imagine 
that part of the platelet is the open field line region, where a characteristic potential 
AV ~ 6.6 X lO^^i^is/P^ volts can be maintained (Ruderman & Sutherland 1995). Here 
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-Bi5 is the surface magnetic field in unit of 10^^ G. Thus, we have 5E ~ 1 x 10^i?i5 esu cm ^. 
Therefore, EqlH] could be re- written as 

5E 

^ 3 X lO^Dl cm s-\ (15) 

where is the radial distance of the emitting region in unit of 10^ cm. By substituting 
the average value -Dg = 0.7, the drift velocity of the emitting region is Vd = 10^ cm s~^, 
which is consistent with our result in FiglTl Furthermore, as illustrated in EqJT^ the drift 
velocity is proportional to the magnitude of the residual electric field, we may also argue 
that when 6E decays as the equilibrium charge distribution in the open field line region is 
being re-established, the drift velocity becomes smaller. Even the emitting region could stop 
drifting and finally stays in the location indicated in panel c of FiglH where it emits the rest 
radiation. If such case applies, it becomes unnecessary for us to propose the sudden change 
of motion direction in FigjSl 



5. Summary and Discussion 

We have investigated the timing properties of SGRs in the persistent state and out- 
burst tails, and ascribed the variety of the pulsed radiation morphologies as the emissions 
coming from the closed field line regions inside the neutron star magnetosphere. For the 
weak relativistic effects in the magnetosphere closer to the stellar surface, the totally dif- 
ferent features of the light curve of SGRs from those of Crab and Vela could be explained. 
Furthermore, by assuming the emitting region drift in the magnetosphere, we are able to 
simulate the pulse profile evolution of SGR 1806-20 during the giant fiare on 27 December 
2004. In addition, when we take the emissions from both magnetic poles into account, we are 
also able to explain the occurrence of the third peak with the phase roughly consistent with 
the observatio n , whic h increases in strength relative to the two major ones (e.g. Figure 2 in 



Palmer et al. I (120051 ) ) at the end of the tails after the outburst. We believe that when the 
emission mechanism becomes clear, the emissivity on the field lines should not be uniform. 
In that case we may be able to predict a stronger third peak. 



Feroci et al. I (120011 ) made simple analytic fits to the fiare light curves of SGR 1900+14 
on 27 August 1998, and concluded that a fireball with contracting surface (rather than a 
cooling surface of fixed area) could provide a reasonable explanation to the decay tail of the 
outburst. They also proposed that the four-peaked profile was produced by several X-ray 
jets tied to the neutron star surface. However, the size of the radiation region in our fitting 
of SGR 1806-20 is unchanging, and the smooth decay in the tail phase of the giant fiare 
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might be due to the coohng of the plasma. It was beheved that the phase stabihty of the 
pulses in light curve is due to the fixed location of the emitting region on the stellar surface. 
But our results show that a small migration could not break the stability in phase. 

We have speculated some possible mechanisms, which cause the radiation region to 
migrate. It seems very clear that the physical motion of magnetic field lines, where the 
charged particles are trapped, must be very slow due to the extremely large drag force by the 
electrons in the core of neutron star. On the other hand, the E x B drift seems more possible. 
However, how this residual electric field survives from the screen of electron/positron pairs 
is not clear. We argue that one possible way is that some platelet is in the open field lines. 
It is unclear if this situation always occurs. 

The 3-100 keV phase-averaged spectrum of the pulsed tail during t he 2004 burst is fitte d 



by a blackbody function at the temperature of 5.1 keV plus a power law (iHurley et al. II2005I ). 
We didn't give the calculated spectrum in this paper, for the radiation mechanism in the 
closed field lines needs further work. However, we argue that the power law component 
results from the synchrotron radiation by the charged particles gyrating along the magnetic 
field lines. We need more information in the higher energy band and the phase-resolved 
spectra to provide more constrains and modification on our geometrical model. 
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Fig. 1. — The closed field lines on the Q-fi plane of layer parameters ai = 5, 10,20. 
inclination angle is a = 30°, and spin period of SGR 1806-20 (P=7.56 s) is applied, 
neutron star is located at the point (0,0). 
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Fig. 2. — Upper: Emission projections on the {(, $)-plane produced by the emitting region 
with azimuthal width of A0^, = 180'^ at layer ai = 15. The sohd hne indicates the hne of 
sight at the viewing angle ( = 50°, and the darker regions correspond to greater intensities. 
Lower: The single sinusoidal pattern of pulse profile corresponding to ( = 50°. SGR 1860-20 
parameters are used in the calculation, and other parameters for fitting are given in Table{T] 
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Fig. 3. — Upper: Emission projections on the {(, $)-plane. The same fitting parameter as 
those in Figl2] except 0=i,c. Lower: The double-peaked pattern of pulse profile corresponding 
to C = 50°. 
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Fig. 4. — Simulated pulsed profiles of SGR 1806-20 in the giant fiare on 27 December 
2004. The emitting region drifts along the azimuthal direction, with the center coordinates 
indicated. The arrows in panel c indicate the positions of the minor peaks of the radiation 
from both magnetic poles. 
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Fig. 5. — Simulated pulsed profiles of SGR 1806-20 in the giant flare on 27 December 2004. 
The arrows indicate the positions of the minor peaks of the radiation from both magnetic 
poles. The emitting region drifts in the direction perpendicular to that in FigJH 
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Fig. 6. — Variation of radial distance of the local position where the arrival photons are 
generated with the pulse phase in FigjH 
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Fig. 7. — The speed of the emitting region drift in 0* direction is calculated by the comparison 
between our numerical results in FigH] and the observation. The speed is given by v^^ {t = 
^^^^) = Ax/ At, where Ax is the displacement of the emitting region centers in two adjacent 
panels in FigJH and At = t2 — ti is the time interv al. At and erro r of t is roughly determined 
from the time scale of Supplementary Figure 1 in lPalmer et al. I (l2005l ). The error of v^^ is 
calculated from the error propagation (j{v^^)/v^^ = a{t)/i. Since the emitting region changes 
the direction of movement or stops moving after t = 160 s, we set Vi^^ {t = 160 s) ~ 0. 
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Table 1. Parameter sets in Fig|2}l5l 



plot label 


typical layer ai 








Figja 
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4.5 


-28 


180 


FigO 
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4.5 


180 


180 


panel a in FigH] 


10 


3.0 


-11 


180 


panel b in FigJH 


10 


3.0 





180 


panel c in FigJH 


10 


3.0 


7.5 


180 


panel d in Figj5] 


20 


6.0 


7.5 


180 



Note. — {9^c, (f>*c) are the magnetic coordinates of the middle 
point of the emitting region, where 9^:c = ai x 9p is the product 
of the layer parameter and the angular radius of the polar cap, 
and (f)^c = 0° is defined as the fl-fj, plane. is the transverse 
extension of the emitting region along the azimuthal direction. 
The magnetic inclination angle is a = 30°, and the viewing angel 
is 50°. 



